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ABSTRACT 

Interactions such as mergers and flybys play a fundamental role in shaping galaxy morphol¬ 
ogy. Using the Horizon Run 4 cosmological iV-body simulation, we studied the frequency and 
type of halo interactions, and their redshift evolution as a function of the environment defined 
by the large-scale density, pair separation, mass ratio, and target halo mass. Most interactions 
happen at large-scale density contrast S ~ 20, regardless of the redshift corresponding to 
groups and relatively dense part of filaments. However, the fraction of interacting target is 
maximum at S ~ 1000. We provide a new empirical fitting form for the interaction rate as a 
function of the halo mass, large-scale density, and redshift. We also report the existence of two 
modes of interactions from the distributions of mass ratio and relative distance, implying two 
different physical origins of the interaction. Satellite targets lose their mass as they proceed 
deeper into the host halo. The relative importance of these two trends strongly depends on the 
large-scale density, target mass, and redshift. 

Key words: Galaxies: haloes, interactions - Cosmology: Large-scale structure of the Uni¬ 
verse, Theory, Dark matter - Methods: numerical 


1 INTRODUCTION 

In a hierarchical ACDM universe, low-mass haloes form first and 
merge together into a more massive one through gravitational at¬ 
traction. Interactions might prevail during the halo formation and 
evolution. Therefore, the star-formation activity and the morphol¬ 
ogy of galaxies are governed by the halo interactions. 

The effects of interactions also depend on the environment. In 
dense environments like clusters, intracluster interactions can sup¬ 
press star formation by multiple high-velocity galaxy interactions 
within the cluster (harassment, Moore et al. 1996) or by tidal strip¬ 
ping (Merritt 1984; Byrd & Valtonen 1990). In addition to grav¬ 
itational interactions, hydrodynamics also plays an important role 
with ram pressure stripping (Gunn & Gott 1972; Bekki 2014; Cen 
2014). Starvation or strangulation strips the hot gas and quenches 
star formation on a larger time-scale (Larson et al. 1980; Cen 2014). 
In particular, minor mergers have been shown to be an efficient 
way to quench star formation, both from theoretical (e.g., Mihos 
& Hemquist 1994) and observational (e.g., Kaviraj 2014a) work. 
On the other hand, mergers have thus been considered as one of 
the most efficient processes to shape galaxies, turning blue, star¬ 
forming, spiral galaxies into red, dead, elliptical ones (Toomre & 
Toomre 1972; Larson & Tinsley 1978). 
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The density-morphology dependence has been studied for 
decades. Rich clusters are dominated by a giant elliptical, and 
the fraction of ellipticals decreases with increasing radial distance 
(Oemler 1974; Dressier 1980; Postman & Geller 1984). On the con¬ 
trary, fields are dominated by blue galaxies (Gomez et al. 2003; 
Goto et al. 2003; Tanaka et al. 2004). 

Park et al. (2008) found — from the SDSS DR4 — that galax¬ 
ies located within the virial radius of their nearest neighbour tend 
to have the same morphology as their neighbour as their separation 
decreases, whereas if their separation is larger than the virial ra¬ 
dius, the fraction of elliptical decreases with increasing separation. 
Subsequent studies at low redshifts (Park et al. 2008; Park & Choi 

2009) , high redshift (0.4 ^ ^ 1.0 Hwang & Park 2009) and low- 

density environment (Yadav & Chen 2014) confirmed these results. 

The Press & Schechter (1974) formalism — and its extensions 
(Bond et al. 1991) — enabled the study of the merger rate from dark 
matter only (e.g., Berrier et al. 2006; Fakhouri & Ma 2008, 2009; 
Hester & Tasitsiomi 2010; Genel et al. 2010) or hydrodynamic 
(Murali et al. 2002; Mailer et al. 2006; L’Huillier et al. 2012; Kavi¬ 
raj et al. 2014; Rodriguez-Gomez et al. 2015; Welker et al. 2015) 
cosmological simulations. However, the merger rate still varies by 
an order of magnitude between different models (Hopkins et al. 

2010) . Environmental studies from SAMs (e.g., Jiang et al. 2014) 
and IV-body simulations (e.g., Fakhouri & Ma 2009) have also 
shown an enhancement of the merger rate in high-density regions. 
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Observationally, the merger rate is usually estimated either 
from close-pair counts (Patton et al. 2002; Keenan et al. 2014; 
Robotham et al. 2014; Lopez-Sanjuan et al. 2015) or by morpho¬ 
logical signatures (Conselice et al. 2003), and converted through a 
merger timescale. However, this merger timescale depends on the 
mass ratio and the gas fractions, which may lead to uncertainties 
in the merger rate (e.g., Lotz et al. 2010a,b). Several studies have 
used visual inspection to study the morphology (Darg et al. 2010; 
Lintott et al. 2011; Kaviraj 2014b). 

JV-body simulations of galaxy encounters have shown that vi¬ 
olent mergers were able to turn spiral galaxies into elliptical ones, 
supported by the random motions of their stars (e.g., Toomre & 
Toomre 1972; Barnes 1988, 1992). The implementation of hydro¬ 
dynamics showed that gravitational torques from the interacting 
pairs could bring gas to the centre of the galaxies (Barnes & Hern- 
quist 1991, 1996), leading to bursts of star formation (e.g., Mihos 
& Hemquist 1996; Tissera et al. 2002; Di Matteo et al. 2007; Hay¬ 
ward et al. 2014). 

However, not all interactions end in mergers. If the relative 
velocity of the interacting galaxies is high enough, they will not 
merge within a dynamical timescales. This class of interaction, re¬ 
ferred to as flyby, has been less studied, although in the recent years 
its interest has been growing (e.g., Moreno 2012; Sinha & Holley - 
Bockelmann 2012; Tonnesen & Cen 2012). At low redshifts on 
high mass scales, flybys becomes as frequent as, or even more fre¬ 
quent than mergers (Sinha & Holley-Bockelmann 2012). Moreno 
(2012) used the Millennium simulation (Springel et al. 2005) to 
analyse interacting pairs that are closer than a fixed comoving dis¬ 
tance dent at some point of their existence. They found that 53% 
of pairs never merge for d a it = 1 /i _1 Mpc, while a lower d C rit 
reduces this fraction. Therefore, difficulties arise in observational 
work, since a non negligible fraction of observed pairs might not 
be gravitationally bound (Sinha & Holley-Bockelmann 2012; Ton¬ 
nesen & Cen 2012). The presence of a hot gaseous halo in late-type 
galaxies is also important for the gas transfer and the star formation 
rate during distant interactions (Hwang & Park 2015). 

A key factor in understanding galaxy evolution, is to under¬ 
stand which of these phenomena (close and distant encounters) 
dominates. Distant interactions are weaker than close ones, but oc¬ 
cur more frequently. It is then important to determine which one is 
the most efficient at shaping galaxy morphology. In this study, we 
used the large Horizon Run 4 A-body simulation to study the en¬ 
vironmental dependency of the interaction rate, regardless of their 
final stage (merger or flybys). 

In § 2 we describes the set of simulations used in this work 
and introduce our definition of interaction and environment. Our 
results about the interaction rate are presented in § 3, and a study 
of the distance and mass ratio is performed in § 4. § 5 discusses our 
results. Finally, the conclusions are drawn in § 6. 


2 SIMULATION AND METHOD 

We present the Horizon Run 4 simulation used here in § 2.1, and 
detail the method we used to construct the halo catalogues in § 2.2. 
Interactions are defined in § 2.4 and § 2.3 describes the density 
calculation. 


2.1 Simulations 

We used the Horizon run 4 simulation (Kim et al. 2015) per¬ 
formed with a memory-efficient version of the GOTPM code (Du- 
binski et al. 2004; Kim et al. 2009, 2011), with a box size L = 
3150 /i _1 Mpc, and A = 6300 s particles. The simulation used sec¬ 
ond order Lagrangian perturbation theory (2LPT) initial conditions 
at Zi = 100 and a WMAP5 cosmology (fib, Om, Ha, h, as, n s ) 
= (0.044, 0.26, 0.74, 0.72, 0.79, 0.96), yielding a particle mass of 
m p ~ 9.02 x 10 9 This starting redshift combined with 

2LPT initial conditions, ensures an accurate mass function and 
power spectrum (L’Huillier et al. 2014). 


2.2 Halo catalogues 

We detected haloes from the simulated particle distributions using 
OPFOF, a MPI-parallel version of the friends-of-friends (FOF) al¬ 
gorithm (Davis et al. 1985), with a standard linking length of 0.2 
mean particle separations, and the gravitationally-stable subhaloes 
with the PSB algorithm (Kim & Park 2006), which are assumed 
to host galaxies. PSB structures identified in an FoF halo are con¬ 
sidered satellite subhaloes, except the most massive one which is 
defined as the main halo hosting the main galaxy. In the following, 
haloes will indistinctly refer to main or subhaloes. 

The virial mass M v of each subhalo is defined as the sum of 
the mass of its member particles. The virial radius R v is then de¬ 
fined by 

3 M v 

47rA c (2)pc(2) 

where the critical density is defined as 

Pc(z) = 3H 2 (z)/(8nG), (2) 

and the critical overdensity at redshift z, A c (z) is well described 
by the empirical Bryan & Norman (1998) fitting formula 

A c (z) = 18 tt 2 + 82a;- 39cc 2 , (3) 




where 


x = 0 ( 2 ) — 1 


_Ha_ 

12 m (l + z) S + 12 a 


(4) 


Note that this definition of the virial radius is in physical (not co¬ 
moving) units. However, in the following, we will only use comov¬ 
ing distances, and we include the Hubble flow Hr in the definition 
of the velocity. 

We defined target haloes (T) as main or subhaloes with Mt > 
ATthrcsh- We chose Mthresh = 5 x 10 11 r'Me, corresponding to 
the mass of 56 particles. This choice yields a mean halo separa¬ 
tion of 5.17/i -1 Mpc, corresponding to SDSS galaxies brighter 
than M r = —19.8 (Choi et al. 2010). For the neighbour haloes 
(N), we limit our catalogue to main or subhaloes with mass Mn > 
qo A/thresh , with qo = 0.4, yielding a minimal number of particles of 
23. The choice of Mthresh and qo are important for the construction 
of our catalogues. A lower Mthresh will result into a larger target 
catalogue, while a lower qo enables to probe larger mass ratios. 
We then chose qo = 0.4 as a good trade-off, yielding Mthresh = 
5 x 10 11 fc _1 Mg. The mean separation of haloes in the neigh¬ 
bour catalogue is 3.97 h~ 1 Mpc at 2 = 0, corresponding to that of 
SDSS galaxies brighter than M r = —18.7 (Choi et al. 2010). Note 
that the target catalogue is a subsample of the neighbour catalogue. 
The minimal number of particles in the PSB subhaloes is set to 
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Lattanzio 1985) of smoothing length h. The latter is searched to 
enclose 20 neighbours, which has been shown to be the smallest 
number yielding an accurate estimate of the density (Park et al. 
2007). The adaptive nature of this technique enables to accurately 
probe the density both in high- and low-density regions (Park et al. 
2007). Appendix A shows the relation between 5 and the cosmic 
web. 

Given the large number of haloes, we adopted a set of oct- 
sibling trees build in rectangular grids over the entire volume to 
efficiently compute the density, similar to Dubinski et al. (2004). 

The local density is defined by the mass of the nearest neigh¬ 
bour and the separation r: 


_ 3Mn _ A c p c 

P" 47rr 3 p 3 


(6) 


and is thus fully described by the ratio of the pair separation to the 
virial radius of the neighbour at a given redshift p = d/R ViIin . 

The mean density is the total mass of the haloes in the Neigh¬ 
bour catalogue 


P=yE Mi - ^ 

This enables us to normalise the density at different redshift by 
considering 


^ = P20 _ i (8) 

P 

We find p(z = 0) = 2.51 x 10 10 h 2 M 0 Mpc -3 . 


Figure 1. Top: Distribution of the large-scale background density at the lo¬ 
cation of the target haloes at different redshifts. Middle: Fraction of target 
haloes that encounter an interaction as a function of the large-scale environ¬ 
ment S. Bottom: Histogram of interactions per density bin. 

20, yielding a minmal mass of 1.8 x 10 11 therefore our 

neighbour catalogue is complete. The target (neighbour) catalogue 
contains 225 406 978 (472 635 985) haloes at 2 = 0. Note that we 
only consider the nearest neighbour, and ignore the effects of a sec¬ 
ond neighbour, which may become important for satellite pairs in 
clusters (Moreno et al. 2013). We discuss this choice in § 5.1 

2.3 Quantification of the environment: large-scale density 
and distance to the nearest neighbour 

The environment may be a key factor to define the interaction rate. 
Both large-scale and local environment have been shown to affect 
galaxy interactions. Muldrew et al. (2012) compared different ways 
to define the environment: using a fixed aperture or the distance to 
the nth nearest neighbour. 

Another way to define and quantify the environment is to com¬ 
pute the eigenvalues of the smooth tidal field (Forero-Romero et al. 
2009) or velocity-shear field (Hoffman et al. 2012). 

We defined the background density as in Park et al. (2007) and 
Park et al. (2008), by 

20 

P 20 = ^2,MiW(ri,h ), (5) 

i=1 

where n and Mi are the distance and mass of the ith closest galaxy 
in the catalogue, and W is the SPH spline kernel (Monaghan & 


2.4 Definitions of interactions 

We define interacting systems as target haloes with mass A/t ^ 
Mthrcsh located within the virial radius of a neighbour more massive 
than goA/x. Such target haloes are considered to be gravitationally 
influenced by their neighbour. 

We define the interaction rate as the fraction of targets that are 
undergoing an interaction over the total number of target haloes. 
This definition slightly differs from some authors, who define the 
merger rate as the number of merger events per descendant (e.g., 
Fakhouri & Ma 2008, 2009, 2010). Indeed, one can notice that 
an interaction may be counted twice if the mass ratio between 
two neighbours is between qo and 1/go, and if their separation is 
smaller than both virial radii (which is likely to happen, since the 
virial radii of these two haloes are comparable). 

We will focus on several quantities: the mass of the halo, its 
background density, the mass ratio q = Mn/Mt and p the pair 
separation divided by the virial radius of the neighbour. We also 
note that, in our definition, the mass ratio q of the interaction is 
by definition larger than go, and the target halo can be the satellite 
(g > 1) or the main halo (g < 1). Our interaction rate corresponds 
to interaction rate per descendant only in the case g > 1. 

Since our parameter space has 5 dimensions, we divide the 
whole target catalogue into subsamples according to the density, 
target mass, and redshift. We the study the two-dimensional (p, q) 
distribution 

d 2 N = f(p,q\S,M T ,z)dpdq, (9) 

where d 2 N is the number of interactions in the interval (p, p + dp), 
(q,q + dq) at fixed redshift z, background density <5, and halo mass 
Mt. 
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(a) z = 0 (b) z = 1 

Figure 2. Distribution function of the target haloes (top left) and interacting target (top right) as a function of the mass and large-scale density at z = 0 (a) 
and 1 (b). Bottom left: interaction rate as a function of the mass and large-scale density. Bottom right: fit to the interaction rate with eq. (10). 
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Figure 3. Dependence of the parameters Ao, M*, and b of equation (10) 
on S and z. 
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3 RESULTS 

3.1 Mass and density dependence of the interaction rate 

In this subsection, we are interested in the frequency of interactions 
that fulfil our definition that a target halo should be located within 
the virial radius of its neighbour. 

The upper panel of Figure 1 shows the probability distribution 
of the large-scale density 8 for all target haloes at different redshifts. 
The distribution range of 8 becomes wider with decreasing redshift, 
reflecting the cosmic evolution of the large-scale structures. The 


Figure 4. Interaction rate (normalised to the final value) as a function of 
the scale factor a for low- (top), intermediate- (middle), and high-mass 
(top) bins, and for low- (stars), intermediate (diamonds), and high-density 
(squares), and the associated fit. 


distribution peaks at 5 of a few units, and the peak height slightly 
diminishes with redshift. 

The middle panel shows the fraction F (5) of target haloes un¬ 
dergoing interactions as a function of the large-scale density 8. In 
low-density regions, and at all redshift, this fraction is very low. 
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z =2.0 



Figure 5. Probability distribution function / as a function of the reduced distance p = d/R vil n and mass ratio q = Af\/AA[- at z = 2. The three columns are 

low-density (<5 < A; i), medium-density (A, i .' A < A, 2 ). and high-density (<5 ^ -*,. 2 ) environments. The three rows are low-mass (Mo < M < Mi), 

intermediate-mass (Mo < M < M 2 ), and massive (Mp M 2 ) target haloes. 


while in high density regions (8 > 100), the fraction monotonously 
rises to reach a maximum at 1 + <5 ~ 500, and then drops again. 
However, at these very-high densities, the number of targets is 
much lower, thus the statistical significance of the drop decreases. 
The interaction rate at a fixed density dramatically increase from 
2 = 5 to 2, and then shows little time evolution after this period. 

Finally, the bottom panel of the figure shows the actual num¬ 
ber of interactions. For all 5, the number of interactions increases 
with decreasing redshift. We find that halo-halo interactions occur 
most frequently in the moderately high-density regions with a very 
narrow density range of 10 < <5 < 30 at all studied redshifts. 
Therefore, even though halo interactions occur in a wide range of 
density environments, the majority of interacting haloes are always 
located in the regions with <5 ~ 20. This corresponds to haloes in 
filaments (see Fig. Al). 

Taking advantage of our unprecedented catalogue, we can 
study the interaction rate as a function of the mass M, background 
density, and redshift separately. 


The upper-left panel of Fig. 2a and b show the distribution 
function of target haloes at z = 0 and 1 as a function of the mass 
and large-scale density. This can be seen as an extension of the 
mass function on a two-dimensional plane. Most targets are located 
at intermediate densities and low masses. At low densities, only 
low-mass haloes are found, while at higher densities, the range of 
masses increases. It should be noted that the mass function has a 
double peak at fixed density. At very high 5, however, the mass 
function is almost flat. The trend of increasing target mass with in¬ 
creasing density is due to host haloes, that follow the large-scale en¬ 
vironment, while the wide range of masses at a given density is due 
to the subhaloes. Even though the definitions of environment are 
different, we find a qualitatively very similar distribution to Zhao 
et al. (2015) and Fakhouri & Ma (2009). The former studied the 
two-dimensional distribution of haloes in the density-Vjnax plane, 
where the maximum velocity V max is a proxy for the halo mass. 
Our correlation between mass and density agree with their results. 
In the latter case, our results are closer to the left panel ( 87 ), that 
also includes the mass of the target halo. 
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z = 1.0 
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Figure 6. Same as Fig. 5, at z = 1. 


The upper-right panel in each subfigure shows the distribution 
of interacting targets. Its shape follows closely that of the targets, 
except for the high-mass end at each density. This is because of our 
definition of interactions: the most massive haloes without a com¬ 
parable neighbour are not considered to be interacting. Therefore, 
the double peaks in the mass function disappears. 

The lower-left panels of Fig. 2a and b show the interaction 
rate, defined as the fraction of target that are interacting for a given 
mass and large-scale density. It is by definition proportional to the 
ratio of the upper-right to the upper-left panels. The rate at a given 
mass increases with increasing 8 , and reaches a plateau at high den¬ 
sities. At fixed density, the interaction rate is constant then drops 
towards large masses. 

At higher redshifts, the distributions look very similar, but ex¬ 
tend to lower masses and densities. 

We proceeded in the following way to fit the interaction rate 
T(M|<5, z). At fixed (8, z ), we used an empirical function defined 
by 

T(M\5,z) = Aoerfc (blog 10 , (10) 


where Aq, AT*, and b are free parameters depending on the large- 
scale density <5 and the redshift. These parameters respectively de¬ 
scribe the total interaction fraction, the characteristic mass at which 
the fraction drops, and the speed of the drop. This function was se¬ 
lected to reflect the plateau-drop shape of T(A/|5). We remove the 
lower and upper density bins that are almost empty of interactions, 
and for which the statistics are very poor. 

We then show the resulting two-dimensional fit of F(A/, 8 ) in 
the lower-right panel of Fig. 2a and b. The fit is overall satisfactory, 
with an underestimation of the rate in the high-<5 end of the range. 

Figure 3 show the dependence of the parameters AT*, Ao, and 
b on the large-scale density at each redshift. We fit them with a 
second-order polynomial in log 10 (l + 8 ), and show it as a dotted 
line. 


3.2 Time evolution of the interaction rate 

In this section, we focus on the time evolution of the interaction 
rate. In order to follow the same biased objects in a statistical sense, 
we proceeded as follows. At a fixed redshift, we divided our sam- 
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z=0.0 
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Figure 7. Same as Fig. 5, z = 0. 


pie of target halo into subsamples of equal number of targets, con¬ 
stant at all redshifts, according to their mass. Then, we divided each 
mass-subsample into three subsamples according to the large-scale 
density of the targets. In practice, at z = 4, we selected all target 
haloes more massive than M 2 (z = 4) = 8.67 x 10 11 /i _1 M 0 , 
corresponding to Nq ~ 4.9 x 10 b . We then divided this sample 
into three subsamples of the same size in terms of the number of 
target halos, namely, 5 < A 2 ,i (z = 4), A 2i i (z = 4) < <5 < 
A 2 , 2 (z = 4), and S > A 2 , 2 (z = 4). At z = 3.1, we found 
M 2 (z = 3.1) and Mi (z = 3.1) such that we have two subsamples 
according to the target mass, with the same number of targets as 
at z = 4. We then divided them into three subsamples according 
to the target large-scale density. At lower redshift, we introduced 
another subsamples at low masses, Mo(z) < M < Mi(z). This 
ensures that the number of targets in each bin of mass and den¬ 
sity is constant, about 1.6 millions, and thus that we are following 
statistically similar haloes across different redshifts. The redshift- 
evolution of Mi and A j } i, (i £ {0,1, 2 },j £ {1, 2}) is shown in 
Fig. B2. 

Figure 4 shows the redshift-evolution of the interaction rate for 


low- (top), intermediate- (middle), and high-masses (bottom panel). 
The low-, intermediate-, and high-density bins are shown as stars, 
diamonds, and squares. For visualisation sake, we normalised the 
rates to be 1 at a = 1. We fit the interaction rate as follow: 

F(a) = Bexp ^ )> (H) 

and the values of A, B, and 7 in each bin are shown in the figure. 
The resulting fit is shown as a dotted line for each subsample. T 
has an inflexion point at a = 1 — A ((7 — l)/ 7 ) 1 ' 7 , which may 
be seen at a ~ 0.3 - 0.4. It is worth noting that this fitting func¬ 
tion is able to reproduce the behaviour of the interaction rate from 
a = 0.2 to 1 and for all range of mass and density considered. Note 
that the fitting function is still positive at a = 0 , so the description 
is not perfect at very low a, but is good enough for the range of 
redshifts of interest. The actual value of the interaction rates can be 
directly read from B. The shapes are very similar, rising then reach¬ 
ing saturation, but the mass and density dependencies can be seen. 
At fixed mass, the slope is steeper in the high-density bin, and the 
rate in the low-density bin is flatter, showing an earlier saturation. 
The final interaction rate significantly (B) grows with the density. 
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while A slowly decreases with increasing density. At fixed density 
bin, the slope becomes shallower and the interaction rate drops to¬ 
ward larger masses. This counter-intuitive point was already seen 
in the right panels of Fig. 2, where the interaction rate drops at large 
masses, and is due to our definition of interactions. Also, the raise 
is somewhat delayed in the largest mass bin. ai/ 2 , the time when 
the interaction rate reaches half of its final value, is larger in this 
bin. 


4 TYPE OF INTERACTIONS 

After studying the effects of the mass and environments on the in¬ 
teraction rate, we study further to investigate the distributions of p 
and q. 

4.1 Distribution of the distance and mass ratio of interactions 

Figure 5, 6, and 7 show the distribution f(p, q\ 8 , AIt,z) at z = 2, 
1, and 0, in the same bins of density and target halo mass as in § 3.2. 
The function is normalised in each bin such that 

//(p,g)dpd g =^" (12) 

J -iVtarget, bin 

over the integration domain, where Aimer, bin and A ta rget, bin are the 
number of interactions and targets in each bin. Note that by con¬ 
struction, Atarget, bin is constant in each panel. The distribution func¬ 
tion is thus normalised to the interaction rate in the bin, as defined 
in § 3 and Fig. 2. There is a clear dependence on the halo mass and 
on the background density. 

One can also see in the lower right panel the existence of two 
main branches: a vertical branch at q < 1 and an oblique branch. 
They correspond to different stages of interactions: 

(i) The vertical branch, with a large range of p for q ~ 1, 

(ii) The oblique branch, where the mass ratio increases as the separa¬ 
tion becomes smaller. This branch can be understood as tidal strip¬ 
ping as the satellites orbits closer and closer to the centre of the 
main halo (e.g., Sales et al. 2007). 

Interestingly, the target mass and background density have a strong 
influence on the relative importance of both branches. 

At fixed redshift and mass bin. increasing the density has two 
main effects. First, the range of q becomes larger at higher densi¬ 
ties. This reflects the fact that massive haloes are located in high- 
density regions. The second effect, that directly follows from the 
first one, is that branch (ii) is shifted towards larger values of q. 

When the redshift and the density bins are fixed, going towards 
higher masses increases the importance of branch (i). Moreover, 
since the target becomes more massive, its virial radius increases, 
and the minimal value of p decreases. 

Finally, at fixed mass and density bin, there is also a redshift 
dependence. Since haloes become more massive, the total range of 
q becomes larger. The importance of branch (ii) fades away with 
decreasing redshift, while branch (i) becomes more and more im¬ 
portant. 

We interpret these two branches as follows: As the interaction 
starts, the satellite halo (less massive one) enter within the host, 
and as it proceeds towards the centre, loses its mass through tidal 
stripping (branch (ii)). After the closest approach, the satellite can 
be anywhere within the neighbour, corresponding to branch (i). In 
the next section, we focus on understanding the physics of these 
two branches. 



Figure 8. Distribution of cos 0 for the different subsamples. 6 is the angle 
of the neighbour’s velocity vector from the position vector from the target. 


4.2 A closer look at the branches 

In order to better understand the physical meaning of these two 
branches, we now focus on the intermediate-density, high-mass 
panel of Fig 7. We consider all interactions with 0.4 < p < 0.6 
and separate then between low- (q < 1, the target is the host of 
its neighbour) and high-q values (q > 3), and now focus on 9, 
the angle between the position and the velocity of the neighbour 
with respect to the target, for the two subsamples. Tangential orbits 
have cos 6 = 0, while receding (approaching) radial orbits have 
cos# = (—)1. 

The left panel of Fig. 8 shows the distribution of cos 6 . For 
a pure random case, the distribution of cos 6 should be flat. Inter¬ 
actions with q < 1 have a flatter distribution of cos 8 , reflecting 
a larger amount of tangential orbits than those with q > 3. It can 
be seen that the distribution of cos 9 for interactions with q < 1 is 
flatter, hence more random. Therefore, the branch at the left of the 
p — q plane corresponds to the haloes whose orbits are randomised 
after they envounter and orbit around their neighbour. On the other 
hand, the interactions with q > 3 have significantly fewer tangen¬ 
tial orbits compared to the random case. This means that the upper 
right branch in the p — q plane is composed of haloes in the course 
of the first encounter in radial orbits. 

The right panel of Fig. 8 shows, for a fixed range of neighbour- 
to-target mass ratio of 1 < q < 3, the distributions for p < 0.4 
(solid line) and p > 0.6 (dotted line). Interaction with p < 0.4, 
which belong to the left branch in the p — q plane, show more con¬ 
sistency with random orbits. The interactions with p > 0.6 show 
relative scarcity of tangential orbits with some excess of receding 
radial orbits over approaching orbits. These observations reinforces 
the conclusion that the left branch is formed by relaxed haloes. 

Figure 9 shows f(p, q\z = 0, Ai ,2 < 5 < A 2,2 ,M > M 2 ), 
the middle panel of the third row of Fig. 7, for cos 9 > 0.5 (left), 
—0.5 < cost? < 0.5 (middle), and cos 9 < —0.5 (right). One 
can see that the left and right panels, corresponding to receding and 
approaching radial orbits respectively, have a prominent upper right 
branch, while the middle row, corresponding to tangential orbits, 
shows fewer interactions in this locus. This again reflects the fact 
that branch (ii) consists of haloes undergoing their first encounter, 
with a rather radial orbit, while branch (i) consists of haloes that 
have already experienced many close encounters. 
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Figure 9. Interaction rate as a function of ( p, q) for interactions with cos 9 > 0.5 (left), —0.5 < cos 9 < 0.5 (middle), and cos 9 < —0.5 (right). 


5 DISCUSSION 


5.1 Choice of the neighbour 

In this section, we study the influence of the choice of the neigh¬ 
bour to define interactions. Park et al. (2008) tested the neighbour 
with the largest tidal energy deposit (LTED) as the most influential 
neighbour, and found essentially no difference. For galaxies in the 
field, this makes very little difference since the closest neighbour 
is more likely to be the only one affecting the target. The situation 
may be different in high-density regions, where two satellite may 
be interacting while the main halo, even more distant, may have 
a stronger influence (Moreno et al. 2013). Park & Hwang (2009) 
showed that while the structural parameters (e.g., concentration) 
depends on the clustercentric distance, the star formation history is 
still a function of the closest neighbour. 

Following Park et al. (2008), we computed the energy de¬ 
posit by each of the 20 neighbours used to compute the large-scale 
density, provided that it fulfils the conditions (Mn > qaM^ and 
r < J?vir,N) using Binney & Tremaine (2008): 


A E oc 


Mia 2 

A v^- + R 2 ^Y' 


(13) 


where a is the rms radius of the target halo, and Av the relative 
velocity. Note that we used the virial radius of the target Rj rather 
than the rms radius, but we expect the results to be unchanged. 
We then defined the neighbour with the largest A E as the most 
influential neighbour. 

Fig. 10 shows the median of the mass ratio q in each bin of 
density. Interactions using the nearest neighbour (NN) are shown 
in red, thick lines, and the LTED case is shown in blue line. The 
error bars show the 15th and 85th percentiles. As expected, the me¬ 
dian of the mass ratio increases at larger densities since it is more 
likely to find a more massive neighbour. Up to 8 ~ 500, the choice 
of the neighbour has no effect. At higher densities, as expected, 
the neighbour with the LTED is not necessarily the NN, but more 
likely to be a more massive halo, which leads to larger mass ratios. 
Since most interactions occur at rather low densities, 8 of a few 
tens, the choice of the neigbhbour has no effect. However, if one 
is interested in interactions within clusters, one should define the 
neigbhour carefully. 



Figure 10. Effect of the definition of the neighbour. Shown are the median 
value and the 15th and 85th percentiles of the mass ratio q for each density 
bin for the default nearest neighbour case (thick blue lines, NN), and for the 
maximal energy deposit (red line, LTED). 


5.2 Comparison with previous work 

Previous theoretical work from simulations (e.g., Fakhouri & Ma 
2008, 2009; Genel et al. 2010) studied the merger rate per time 
unit using merger trees, and typically found a merger rate of 
dfV m erg/dt oc (1 + «)“, while Fakhouri & Ma (2009) found an 
environmental factor (1 + (M/M 12 ) 7 , where 5 7 is the local 

density computed over a sphere of radius 7 h - 1 Mpc. However, our 
approach differs from previous theoretical studies in that we do not 
use merger trees to define interactions, but rather look at the dis¬ 
tance to the closest neighbour at a given time. This makes direct 
comparison with previous work difficult. Indeed, in our approach, 
the same interaction may be counted at different redshifts, while in 
other studies it will only be counted when the merger is complete. 

Our approach is by construction closer to observations in that 
we look at the instantaneous interaction fraction, namely, the frac¬ 
tion of target that are undergoing an interaction at a given redshift, 
corresponding to f m in (Lotz et al. 2008). However, our condition 
that the target should lie within the virial radius of its neighbour, 
increases our interaction numbers especially at low-redshifts where 
the physical virial radius of the studied mass range is typically 0.5 
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to 1 h 1 Mpc, compared to close-pair counts that limit the separa¬ 
tions to a few 10 /i _ 1 kpc. 

The main advantage of our approach, besides its simplicity, is 
that it does not require observations to use a merger timescale, that 
can be a source of uncertainty. Instead, the observed merger frac¬ 
tion can be directly compared to our fitting formulae. However, one 
should keep in mind that, since we consider distant interactions, our 
pair separation can be larger than the typical 30 /i _1 kpc separation 
used to calculate merger rates. 

Observationally, most studies define the interaction radius, 
i.e., the largest radius to define an interacting pair, as fixed in phys¬ 
ical scales. Lopez-Sanjuan et al. (2015) compiled several study of 
close pairs for 0 < 2 < 1 . 2 , with a minimal projected separa¬ 
tion of 100 /t _1 kpc (physical), and fitted an interaction fraction of 
f = /o(l + z) m . They found that /o decreases and m increases 
with decreasing magnitude cut. While our fitting formula are dif¬ 
ferent, fo and m can be related to our B and 7 parameters. How¬ 
ever, their behaviours are quite different, since B is constant with 
mass, except in the largest density bin. The main difference with 
our approach is the use of a fixed projected radius for the defi¬ 
nition of interacting pairs, while our interaction distance depends 
on the neighbour. At fixed mass, and in proper units, the virial ra¬ 
dius increases with decreasing redshift (since R vn and A c decrease, 
c.f. equations 1 and 3). Moreover, haloes also accrete mass, which 
also makes the virial radius grow. These considerations results in a 
larger interaction fraction at lower redshift in our definition. 


6 SUMMARY 

In this study, we used the massive Horizon run 4 cosmological 
IV-body simulation to investigate halo interactions, by considering 
three aspects of the environment, namely, internal (target mass), lo¬ 
cal (distance to the nearest neighbour), and large-scale (background 
density). 

This simulation has the largest size among all simulations 
that have been used to study galaxy interactions, which enables to 
study the environment dependence and very good statistics, while 
its resolution is high enough to study haloes with masses down 
to 5 x 10 11 h-'Mo. We defined target and neighbour haloes as 
haloes more massive than 5 x 10 11 and 2 x 10 11 and we 

define a target of mass Mt to be interacting if it is located within 
the virial radius of a neighbour more massive than OAMt- 
Our findings are summarised as follow: 

• Interactions preferentially occur at intermediate large-scale 

background densities, 8 ~ 20 , regardless of the redshift, and for 
targets more massive than 5 x 10 11 The interaction rate, 

i.e., the fraction of targets undergoing an interaction, reaches its 
maximum at larger values: 8 ~ 1000 . 

• We propose a new functional fit (eq. 10) for the interaction 
rate as a function of the target mass and large-scale background 
mass density at fixed redshift, and give the redshift-evolution of the 
fit parameters. 

• The interaction rate as a function of time has a universal shape, 
rising and saturating. We propose a fitting formula for the time evo¬ 
lution of the interaction rate (eq. 11). This function strongly de¬ 
pends on the large-scale density and target mass. Larger density 
yield a steeper slope and larger final interaction rate, and massive 
haloes saturate later. 

• We reported the existence of two branches in the two- 
dimensional distribution f(p, q\8 , Mt, z), namely: an oblique 


branch where the distance to the neighbour decreases as the mass 
ratio increases, reflecting the loss of the satellite mass through tidal 
stripping as it orbits within the main halo. A second branch at lower 
p for the same q corresponds to interactions that already passed the 
closest approach and whose orbits became randomised. 

• The relative importance of the two trends strongly depend on 
the large-scale density, target mass, and redshift. 

In a paper in preparation, we will study the effects of inter¬ 
actions on parameters such as the alignment of the spin and major 
axes of the interacting haloes, and their orbits. In a next step, we 
will add hydrodynamics to study the galaxy (rather than the halo) 
interaction rate, which can be significantly different (Hopkins et al. 
2010 ), and to probe the effects of gas physics on the morphological 
transformation of galaxies. 

We argue that observers may use our definitions to measure 
the interaction fraction. In a next paper, we plan to apply our 
method to publicly available observational data in order to com¬ 
pare the evolution of interaction fraction between observations and 
our theoretical results. 
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APPENDIX A: 5 AND THE COSMIC WEB 

It is interesting to understand the relation between the overdensity 
parameter S and the cosmic web, namely, where haloes with a given 
S are located within the cosmic web. Figure Al shows a 150 x 
150 x 20/i~ 1 Mpc slice at z = 0 (top) and 1 (bottom). High-J 
haloes are shown in large, red circles, intermediate-^ (10 < 1 + 
S < 30) in medium, cyan circles, while other haloes are shown 
in small, blue circles. Haloes are overplotted on top of the density 
map computed with the Triangular-Shape Clouds, with a pixel size 
of 0.5 /i _1 Mpc, the same as used for the PM calculation. At z = 0, 
haloes with 5 > 100(red) can be found in nodes, while intermediate 
(10 < S < 30) are found in filaments. On the other hand, at z = 1 
the overall density fluctuation amplitude is lower, and the haloes 
with 5 > 100 are located only in the rarer highest density nodes. 
Since we vary the density threshold to cope with the evolution of 
the growth factor, the type of environment corresponding to each 
density bin is supposed to be independent of redshift. We conclude 
that the overdensity 1 + <5 is a good proxy for the environment. 
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Figure Bl. Same as Fig. 7, with z msa = 2. 


APPENDIX B: EFFECTS OF M 0 , M u AND M 2 

Fig. Bl is the same as Fig. 7, but using 2 max = 2 instead of 4. The 
larger number of targets available from z = 2 enables us to keep 
better statistics, but prevents from doing a redshift-evolution study 
as in § 3.2. It results in a much higher number of targets in each 
panel, about 20 millions, from z = 2 to 0. 

The value of Mo, Mi, and M 2 are respectively 6.86 x 10 11 , 
1.1 x 10 12 , and 2.14 x 10 12 compared to 7.68 x 10 12 , 

1.10 x 10 13 , and 1.99 x 10 13 /t _1 Mg for the z m nx = 4 (shown in 
Fig. B2) case, which implies that all targets in Fig. 7 are in the lower 
panel of Fig. B 1, resulting in a more important oblique branch. 
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Figure Al. Distribution of haloes with different overdensities in the cosmic 
web at z = 0 (top) and 1 (bottom). The grey background is the density 
field computed over 0.5 fi. -1 Mpc pixels with a TSC method. The larger 
values (large red circles, 1 + <5 > 100) are located in the nodes, while 
intermediate valuers (medium, cyan circles 10 < 1 + <5 < 30) are located 
in the filaments. 
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Figure B2. Redshift evolution of Mj (top), and A ij (bottom), where j 
stands for the mass bin and i for the density bin, for the z max = 4 case in 
§3.2. 


















